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Abstract 

The energy dissipation in a gas of structured objects, e.g. molecules, is considered 
in density matrix formalism. It is shown that the macroscopic irreversibility of the 
kinetic processes can be considered as a consequence of the microscopic operator or- 
dering. Our approach is free of any special assumptions on the space-time geometry, 
except for the general causality assumptions, so it can be applied to a wide variety of 
processes, from the cosmological processes at Big Bang stage till the energy dissipation 
in molecular gases. 

1 Introduction 

Kinetic equations describe the evolution of the distribution function of matter in phase space 
at the presence of particle coUisions. They lay the foundation of many cosmological models, 
describing the formation of present barionic matter from initial plasma. The problem with 
kinetic approach is that at the Plank times, r ~ 10~^'^sec, the era of quantum gravity, the 
Riemann space-time itself did not exist, so neither the Boltzmann kinetic equation of the 
form 

df{X,P,t) , ^ df{X,P,t) , P df{X,P,t) 

di + ^ dP + m dX = 

nor its relativistic generalisation is valid for that stage. 

The situation with the quantum mechanical description of the Early Universe is a little 
better. In quantum mechanics the evolution of quantum system is described as the evolution 
of the density operator p obeying the von Neumann equation 

zhp=[H,p], (2) 
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where H is the Hamiltonian of the system. For the evolution equation ([2]) it does not mean 
whether the spacetime is continuous, differentiable, Archimedian. What is significant is the 
evolution: the existence of time t and the dependence of the state on t just means the 
evolution; if there is no evolution there is no change in quantum state. 

It is important for the kinetic approach that the evolution equation for the classical 
distribution function in phase space f{X,P,t) should be obtained as the classical limit 
^ — > of a more fundamental quantum equation ([2]) at the assumption of the existence of 
classical trajectory X = X{t) for each particle. This is called the macroscopic limit. 

In a purely quantum case, i.e. when h cannot be treated as a small parameter, it is 
impossible to assign any trajectory X(t) to a given particle. This is the case for dense 
plasma, lasers, microelectronics etc. This is also the case for quantum gravity. 

The present paper considers the effect of general causality assumptions on the evolution 
of the density operator and the kinetic equation, resulting from this evolution in macroscopic 
limit. The general causality assumptions, discussed in pjE], are required for the construction 
of quantum field theory on a category space, which is not a manifold, i.e. when continuity 
and differentiability are not provided. This is the case for the quantum gravity era of the 
Early Universe, and hence before the light-cone causality of Minkovskian space have been set 
by Big Bang a more general causality assumptions should hold. In the limit of classical non 
relativistic gas the general theory drives us back to known results for the Boltzmann equation 
and the Klimontovich method for molecular gases, but ultimately leads to irreversibility if 
the internal degrees of freedom are involved. 

The remainder of this paper is organized as follows. In Section we remind the links 
between quantum description of hierarchic systems and causality. In Section a semi- 
quantum kinetic equation for a gas of structured particles is constructed. [7] presents 
a toy model of a system of hierarchic quantum oscillators and imposes the operator ordering 
on this system. In Conclusion we summarize the basic ideas of our approach. 



2 Density operator for hierarchic structures 

Strictly speaking the von Neumann equation ([2]) holds only for closed systems - but the only 
system which is absolutely closed is the Universe as a whole. Therefore there should be some 
methods to treat the system of our interest as approximately closed, when its interaction with 
the rest of the Universe can be neglected, or treated perturbatively. In this case the state 
vector of the whole Universe can be casted in a form 

|^) = 5Zc,„|0,)|^„), (3) 

ia 

where are the states of the system, 16*0,) are the states of its environment, i.e. the rest 
part of the Universe. To get a tractable model the states of the environment are assumed to 
be mutually orthogonal 

(OalOp) = 5a0. (4) 
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In view of this assumption the matrix elements of the density operator of the whole Universe 



P=|V^)(V^| 

in the basis of the states of system are 

a 

This means the probabilities of different states of a quantum system i are obtained by 
averaging over all states of its environment. 

The mean value of arbitrary physical observable A for such system is equal to the trace 
of the product Ap with respect to the states of the system 0: 

{A) = TtAp = Y,A,P%. (6) 

Since it is impossible to account for all degrees of freedom of the environment, what 
is done in practice is different from the summation <^ with orthonormality condition (jl]). 
Namely, considering a diluted gas of diatomic molecules in coordinate representation [Sj H], 
one uses the Wigner function 

p{R, P, r, r', t) = -^^ I p,{R-^^,R+ ^7, r[, r^, t)e-^^rf7, (7) 

where R is the center of mass coordinate of the molecule, P is its total momentum; r = 
Vi — r2 is the internal degree of freedom - the relative displacement of atoms in the molecule; 
Pi{f i,f 2, 1^1,^2, t) is the density matrix of diatomic molecule in coordinate representation. 

In the language of state vectors the Wigner function formalism corresponds to the change 
of basis Q for the density operator ([S]) to a hierarchic one 

p=\^){^\, i^)= 5^ Q„.„_,i0?ji0r-')' (8) 

^a — 1 

where (a — 1) denotes the next to the studied system (0°) hierarchy level. (Say, in the case 
of diatomic gas |0°~^) denotes the state vector of molecule, with |0") being the state of atom 
in this molecule.) The density matrix for the system (0°) can be casted in the form [5]: 

/5° = Tr,_i|^)(7A|, (9) 

where {ip) is given by ([8]). 

The hierarchic representation ([8]) poses the problem of ordering the operators acting at 
different hierarchy levels. The operator ordering rule "the coarse acts the first" was suggested 
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Figure 1: Creation of the discrete Universe from vacuum as a cascade process 

by the author [2J. This corresponds to general cosmological idea of cascade process of the 
Universe origin from vacuum: 

\Uo) = al\0), \UoUo^) = alal\0) = {\Uo), |f/o.)}, .... (10) 

At present stage of the Universe evolution the hierarchy level and the number of constituents 
are very large iV ~ 2^ is expected to be of order of the Dirac-Eddington number ~ 10^°, 
or so, and we observe continuous Universe |6l [7] . 

This implies two types of causahty: (i) the vertical causality a C 6 (a is a part of 6); (ii) 
the horizontal causality a ^ b {a precedes 6); in Fig. [T]t/oi C Uq, Uqi -< Uoj. Consequently, 
two types of operator ordering are required to meet two types of partial order. The approach 
with two types of causality relations is a generalization of common light-cone causality of 
the Minkovsky space with T-ordering relation for the operators on more general spaces. 
It perfectly meets the needs of quantum gravity level cosmological models, when the the 
manifold structure is not provided [HI [6]. The details of this topological approach, known as 
region causality, can be found in [1]. 

What region causality says is quite natural not only for quantum gravity, but also for 
molecular physics: no experiment can detect a causal effect in a point - what is really 
observed is an effect in a region, the size of which is restricted, at best, by the Heisenberg 
uncertainty relation ApAx > |. So, in accordance to [1] we adopt two types of causality: 

a (1 b a is a part of b 
a b b can see a. 

The latter corresponds to the light-cone causality; the former, supplied with the operator 
ordering "the coarse acts first", just stands for the fact that it is impossible to change the 
state of the part without changing the state of the whole. 
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It is important for kinetic theory that the Bogolubov causahty relations for the scattering 
matrix [5] 
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5g{x) 



0, x^y, (11) 



originally formulated for the light-cone causality - x < y means Xq < yo and x ^ y means 
the interval between the events x and y is space-like interval - without significant changes 
are generalized to the vertical causality as well. This can be easily seen by introducing p-adic 
coordinate x G Zp on the branching tree. The research in quantum gravity have stimulated 
the extension of the ^-matrix formalism to non-Archimedian numeric fields, namely to the 
field of p-adic numbers [lOj. Hence the Bogolubov causality relation flTTl) can be given for 
the quantum fields depending on p-adic arguments, and therefore the fields defined on the 
vertexes of p-adic tree can be ordered in analogy to the T-ordering causahty. 

In p-adic metric the partial order is given by p-adic metric | ■ \p: 

X < y ii \x\p < \y\p - vertical ordering; 

if \x\p = \y\p then x and y are ordered by the first nonzero coefficients - this is horizontal 
ordering in the ring of natural numbers pUl [TT] . So, ii x,y G Zp in f llip the Bogolubov 
causality relation is defined on the branching tree. A toy-model of such field theory is 
presented in |T2]. 

In topology the causality of regions is axiomatising by means of partial order relations 
in a form of causal sets [S], US]. The causality at quantum gravity level - we would say 
before the Big Bang - will be the region causality with two causality relations A d B and 
A -< B , that is called a causal site. The axiomatics of causal site given in [1] is presented 
in Appendix. 

In a (pseudo-) Euclidean space the structure of the causal sites implies causal paths and 
the geodesic coinciding with light-cone causality. However at quantum gravity level, when 
the Archimedean axiom does not hold, the vertical causality relation C is expected to play 
the important role. We do not know the concrete scheme of the present Universe formation 
from initial object, but without loss of generality we can assume that at initial time t = 
there was only one initial object Uq. At the next instant of time t it branched into p > 2 
parts, each of those continued further branching in a tree-like matter - therefore a hierarchic 
structure have been formed This hierarchic structure endowed with certain system of 
relations between its elements has formed our space-time. 

The axiomatics of causal site is rather general and should be applied not only at quantum 
gravity level, but also for quantum systems of present life, such as molecular gases. To see 
this, we have to cast the density operator in a hierarchic basis ([8]) . Let us consider a 
hierarchic system O consisting of two parts I and J, each of those consists of its own two 
parts, ii^i') and respectively, see Fig. [21 For definiteness, let us consider the particle 

"i" of the finest hierarchy level of the system above. The mean value of an observable A 
related to i can be evaluated in two different ways: 
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Figure 2: Hierarchic binary system. Two means of calculating partial density matrices of 
the sub-parts can be applied 

1. with the density matrix averaged over all states of i': 

Pi = ^(ApV)'^ 

i' 

2. with the density matrix averaged over the states of / - the "parent" of i: 

I 

The former case is a standard way of the partial density matrix for a multi-particle system 
[3], the latter is presented in [5]. In the latter case the state vectors are represented in a 
tree-like form 

|<|.,) = {|/),|/^),...}, p=|$,)($,|, (12) 

with index O being dropped as omnipresent. 

To provide the equivalence of these two approaches the effect of the other particles J ^ I 
on the parts of /, i.e. on should exist only via change the state of /. This the effect 
of the whole to its parts. In thermodjTiamics it corresponds to adiabatic insulation of the 
system /, when the energy can be transferred to the system as a whole, but not to its internal 
degrees of freedom. 

3 Quantum kinetic equation for the gas of oscillators 

Let us consider a simple example. A homogeneous gas of molecules with the distribution 
function f{P). The molecules can transfer part of their kinetic energy to their internal 
degrees of freedom, represented by quantum harmonic oscillators. The molecules do not 
interact beyond the interaction zone \xi — 332! > 2ro, with tq being the typical radius of the 
molecule. This approximation is valid since the energy of oscillatory degrees of freedom is 
much bigger than that of rotational degrees of freedom, but much less than the energy of 
electron transitions: 

Eel ^ Eosc ^ Erot- (13) 
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So, in a certain range of energies, corresponding to infrared radiation energy, the kinetic 
energy of molecules is transferred into oscillatory degrees of freedom by means of binary 
collisions between molecules 

The energy levels of quantum harmonic oscillator, representing the oscillatory degrees of 
freedom of the molecules, are 

En = huj(^n + ^ , n = 0,l,... (14) 

The state of each molecule in this semi-quantum approximation is given by three variables 
[X,P,n) - position, momentum and excitation number. The rotations are ignored in this 
approximation f[T^ . The gas is assumed to be homogeneous (to be in a constant force 
field), therefore the total momentum is conserved in each binary collision, the energy can 
be transferred to the internal degrees of freedom only by quanta AE = hu, but not in 
a continuous manner. To make the model consistent we assume the oscillator degrees of 
freedom to be in thermodynamic equilibrium with radiation at temperature T. 

Taking into account that in quantum description ([2]) the density matrix can be casted in 
terms of the Wigner functions, we can derive the kinetic equation for the probability density 
functions /„(X,P,t) with the initial condition fn^Q{X,P,0) = 0. The causality condition 
"the coarse acts first" will manifest itself in the restriction that the energy transfer from the 
macroscopic {X, P) degrees of freedom to the microscopic n is allowed, while inverse is not: 
the excess of energy is relaxed into radiation only. 

Let us denote the momenta of two colliding molecules (in centre mass system) before 
the collision as P and Pi, and P' and after collision, respectively. The momentum and 
energy conservation laws for binary collision imply 

P + Pi = P' + Pi 

As a simplest model the assume the probability of transfer of quanta of energy into 
internal degrees of freedom is proportional to the number of means g{N) this energy can be 
distributed between two oscillators: 

E^^sWfl^O^^, ,(/V) = iV + l. AE. = n.N. (16) 

E«r„'s(A') 



The sums in (1161) should be bounded by some reasonable integer Nmax > 0, which guaranties 
the condition (fT3l) holds. 

The collision integral for our model, i. e. the change of particle density per unit of phase 
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space volume, is given by 



Icol — 


j dTidT'dT[w'{f'f[P{-AQ) 


-fhP{AQ)) 


dT = 


dXdP 

{2'Khf ' 




w' = 


w'{P\ P[ ^ P, Pi) ~ 6{P' + 


P[-P- Pi) 


f ^ 


fn'iX', P', t), 




AQ = 


p2 p2 p/2 p,2 

2m 2m 2m 2m 





(17) 

The detailed balance principle does not hold for the kinetic equation 

dfnjx^p^t) dux,p,t) p dux,p,t) ^ f(xpf)-T \f] ns^ 

I 

for P{AQ) ^ P{-AQ). The coefficients Kni account for the interaction of the internal 
degrees of freedom with radiation. 

If the internal degrees of freedom, i.e. the oscillators (fTi|) . are in thermodynamic equi- 
librium with radiation at temperature T, the ratios of the energy level populations are given 
by 

"'■-expf-^V (19) 



where Uj = J fj{X,P,t)dr, - is the Boltzmann constant; for definiteness we assume 
Ej > Ek. In thermodynamic equilibrium the transition rate from j to k is equal to that 
from k to j. Since the population of each energy level is kept stationary 



dn^ 



kjk'iT'k kj^jTij — 0, kkj — kjk, (20) 



where the first term stands for the excitations k ^ j and the second - for relaxations j k. 

For molecular oscillators the typical frequencies are infrared Ujk ~ 10^'^sec^^. In dipole 
approximation (see e.g. [T5] ) 

_ 4cj|fcrf2^ 1 

exp i^j^j - 1 

where djk is the matrix element of dipole moment operator for the j ^ k transition. For 
molecular oscillators the equilibrium transition rate (12T|) has the order of magnitude lO^sec"^. 

For our case of the relaxation of the kinetic energy of macroscopic degrees of freedom 
into the microscopic oscillators the spontaneous and stimulated emission is of paramount 
importance. Spontaneous emission rate A^j is the probability of downward transition j —>■ k, 
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extra to the equilibrium one k^j fl20l) . In thermodynamic equihbrium the total numbers of 
upward and downward transitions are equal 



kjkUk = {hj + Akj)nj (22) 

from where, with the symmetry relation kjk = kkj and the equilibrium populations f|T9l) . the 
spontaneous emission rate is found to be 

A,. = (23) 

For the low temperatures, fiujjk ^ ksT, the spontaneous emission dominates over the equi- 
librium transitions 

Ajk ^ kjk- 

For typical molecular frequencies the order of magnitude of the spontaneous emission rate 
is as high as Ajk ~ 40sec~"'^ [.15j. 

If the time before collisions is much less than spontaneous emission time 

1 

which is often so in molecular collisions, than the macroscopic collision plays the role of the 
measuring device with respect to quantum oscillator degrees of freedom, i.e. collisions force 
the spontaneous emission - that is they result in stimulated emission. 

The direct, i.e. the photon-less, energy transfer of the excitation energies of colliding 
molecules into their kinetic energy seems to be impossible. If the molecules are in excited 
state j and during the collision they relax to the lower energy states k the energy Ej — Ek 
of the molecules is released in the form of photons with frequencies 

Ej — Ek 

^jk — 



h 

and the momentum 

^jk Ej — Ek 
Pjk = — = — , 

c c 

see Fig. [31 Due to the conservation of total momentum, the momentum of each molecule is 
changed by —Pjk, and its energy is increased by 

AE^^ 
2m 

(in non-relativistic approximation, since momentum transfer is small). 

Suppose now, the photons are virtual, rather than real. If so, the photon emitted by 
molecule 1 is absorbed by the molecule 2, and drives it back to the excited state j; the 
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Figure 3: Simultaneous energy relaxation by two molecules in a binary collision 

photon emitted by the molecule 2 is absorbed by the molecule 1 and drives it back to the 
excited state j. In such a way we have increased the kinetic energy of colliding molecules 
without changing their internal states. This is a kind of perpetum mobile and is therefore 
forbidden. In the language of hierarchic state vectors ([H]) such process of energy transfer 
from 10") to is prohibited by the operator ordering rule "the coarse (a — 1) acts first". 

Similarly, the molecules in Fig. [3] may be two mesons in excited states. In that case the 
quark confinement prohibits the transfer of internal excitation energy into the meson kinetic 
energy. 

4 Oscillator toy model 

In previous section we considered a gas of molecules where the kinetic energy of colliding 
particles is redistributed into internal excitation energy of particles by means of inelastic 
collisions. The particle motion in phase space was described classically, but the internal 
degrees of freedom were treated as quantum oscillators. It is possible to construct a toy 
hierarchic model treating where all degrees of freedom are quantized. 

Let us consider a system of equally spaced particles of mass mo located at x = Sq, s°, s^, . . ., 
such that each (2z)-th particle is connected to its neighbour {2i + l)-particle by a harmonic 
potential of rigidity ko. The hierarchic structure is constructed in a way that two "atomic" 
particles. {2i,2i + 1), form the i-th molecule of mass mi = 2mo; the molecules, in their 
turn interact to each other in the same hierarchic way as atoms, but with different rigidity 
constant ki, see Fig. HJ The Hamiltonian of the hierarchic two- level system shown in Fig. H] 
has the form 




(24) 



where we have introduced the hierarchic centre mass coordinates 




J 

'2i+l 



(25) 



and the displacements 




2i 



s: 



'2i+li 



(26) 
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Figure 4: Toy model system of hierarchic oscillators 



so that the inverse transform is 



si - 



di 



The straightforward algebra gives 

or, after simplification, 
H — mo 



i=0 



1=0 



i=Q i=0 



(27) 



Since Sq — because of the zero total momentum, the final equation for the Hamiltonian is 
a sum of two independent oscillators, belonging to "0" and "1" scales, respectively 



H — mo 



i=0 i=0 



(28) 



For each level of the hierarchy it is possible to substitute creation and annihilation operators 
, a and the number of excitation. In this way the scattering of two molecules described by 
(ig coordinate will result in relaxation of the excitation and redistribution of this energy into 
the modes dQ,d\, so that 

HuiAnl = huJi{Anl + An}). (29) 

In classical toy models, like the Toda chains of oscillators there is no such independence 
of oscillations of different scales as in our model considered above. If the particles are 
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allowed to interact with all the neighbours, including that in the next block, the interaction 
Hamiltonian will have an extra term. Say, for the system shown in Fig. H] this will be the 
term 



rr _ h S-^, _ 
^ib — cy _^l''2j+l ■='2: 



2 

For the first two blocks this gives the contribution 

A^.. = |(^?-^^)^ = |(rfg-^)', (30) 

which produces the inter-scale energy transfer. 

The Hamiltonian (|28l) does not contain any cross-scale interaction terms. Therefore the 
equation (l29il is just a formal equation for energy conservation without any specification of 
the processes of energy exchange. If the cross-scale interaction fl30l) is introduced, the free 
oscillators can be casted in terms of creation and annihilation operators a^, a and ^ A 



where /io = Tno/2. The operators a and are related to the normal modes = a/A^i/ttiq, uj = 
\/ko/fio coordinate and momenta 



,1 h tti -\- al 1 /T ctj — al 

d^ = \ ^ , = ^hllQUJ— . 

y (JLqUJ 2 ll 

The inter-scale interaction Hamiltonian (l30ll in second quantization formalism has the form 



knh 



ib 



A + 1 1 + aj 



2 L(A;imo)i/4 2 2 (Vo)^/^ .^^ 



TTlE^l' (32) 



If we assume the scale ordering - an analog of T-ordering, - then only the terms of the form 
a] A will contribute to cross-scale interaction. 



5 Conclusion 

In classical theory of molecular gases the kinetic energy of molecules can be transformed into 
the heat energy of their internal degrees of freedom and vice versa: the energy accumulated 
in the internal degrees of freedom can be relaxed into kinetic energy of scattering particles, 
see FigJSl In purely quantum case, when momentum and position become operators, the 
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Figure 5: Classical accumulation of mechanic energy: the potential energy of contracted 
springs of two oscillators can be relaxed in the form of kinetic energy of scattered oscillators 

direct energy transfer from internal (microscopic) to macroscopic degrees of freedom seems 
to be impossible. 

The mathematical reason for that is the need to order the operators related the part and 
the whole (the micro and the macro) degrees of freedom. Apart from the requirement of 
finiteness, which is not discussed in this paper, this need stems from the fact that the state of 
the part cannot be affected unless the state of the whole is affected (measured) first. This is 
axiomatized by the operator ordering rule "the coarse acts first" [2J. The set of axioms with 
two partial ordering operations [1], aimed for the development of quantum gravity theory, 
may or may not be self-evident at common scales / S> Ipi- However, there exists a simple 
example illustrating the impossibility of the process, shown in Fig. [5l for a pair of quantum 
oscillators. 

Indeed, in classical system the potential energy is accumulated into the springs by de- 
creasing the size of oscillators, the relaxation of potential energy into kinetic energy of 
scattered particles goes through the expansion - the increase of oscillator sizes. For the 
quantum oscillator in the n-th excited state the mean squared size is 

Thus the relaxation from high to low energy level results in "contraction" of oscillator, rather 
than in expansion. So, no scattering takes place unless an extra momentum is injected. 

In the final end the axiomatizing of operator ordering used above provides a reformulation 
of the second thermodynamics's law on the microscopic level. For a microscopic system, 
described classically, the heating of internal degrees of freedom cannot result in macroscopic 
movement if it is accompanied by the decrease of entropy. In quantum case the decrease of 
entropy can be prevented by operator ordering. 
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A Axiomatic structure of causal sites 

A causal site is a set of "regions" with two binary relations denoted C and -< satisfying the 
axioms: 

1. VA,S,C 

{a)ACBABCC^ACC, 

(b) AC A, 

(c) A C B A B C A ^ A = B. 

2. The partial order C has a minimum element (p, called an empty region. 

3. The partial order C has unions: \/A, B3A U so that 

(a) AC AU B A B C AU B 

(h) A C C A B C C ^ AU B C C. 

4. The partial order -< induces a strict partial order on the nonempty regions: 
{a)A^BAB^C^A^C, 

(h) Ay^A 

5. WA,B,C 

(a) A C B AB ^ A^C, 
{h)A CBAC^B^C^A, 
{c)A ^CaB^C^AuB^C. 

6. VA, B3Ba, cutting of A by B, so that 

(a) Ba-< AABaC B, 

{h) D -< A A D C B ^ D C Ba. 

7. If ^4 -< C are nonempty regions, and 3D : A ~< D ~< C, then 3B complete with respect to 
A^C 

Dei.: li A ^ B ^ C, then B is said to be complete with respect to a causal pair A ^ C, 
if any causal path from A to C can be refined to a causal path from A to C, on of whose 
members is contained in B. 

Def.: A causal path P is a sequence of nonempty regions Ai ~< A2 -<...-< An. 
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